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Abstract 

The lowest order contribution of the amplitude of the Dirac - Coulomb 
scattering in de Sitter spacetime is calculated assuming that the initial 
and final states of the Dirac field are described by exact solutions of the 
free Dirac equation on de Sitter spacetime with a given momentum and 
helicity. One studies the difficulties that arises when one passes from the 
amplitude to cross section. 



1 Introduction 

It is known that many difficult problems arises when one tries to combine the 
quantum theory of fields with general relativity. One of the main problems is 
related to the physical interpretation of the one-particle quantum modes that 
may indicate how to quantize the fields. The curved speacetime have specific 
symmetries which, in general, differ from that of the Minkowski spacetime. For 
this reason, here the symmetries generating conserved quantities have to be 
treated with specific methods. In other respects, it is known that the form 
of the fields equations and implicitly their solutions on curved spacetime are 
strongly dependent on the tetrad gauge and local chart in which one works. It 
is important to point out that in the latest years important steps in developing 
a quantum field theory on de Sitter spacetime was made by finding analytical 
solutions for the Dirac equation in moving or static local charts suitable for 
separation of variables [2] , [8] , [9] , [10] . 

In the present the majority of investigations dedicated to Q.E.D on curved 
spacetime do not take in considerations scattering processes. Our aim in this 
paper is to calculate one scatter amplitude for Dirac field in curved spacetime 
and to analyze the physical consequences that emerge from this calculation. 
Differences with respect to the Minkowski case are that the modulus of the 
momentum is not conserved and a tendency for helicity conservation. Also in 
the limit of a vanishing expansion rate of the space the Minkowski amplitude 
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will not be recovered. Finally, we obtain that the cross section is a sum of two 
contributions, that corresponds to linear and nonlinear amplitude of scattering. 

In section 2 we derive the formula for the scattering amplitude. In section 
3 we discuss some of proprieties of scatter amplitude and we will examine the 
limit cases of a small/large expansion rate of the space compared with the 
mass of particle, in section 4 we calculate the cross section. Our conclusions 
are summarized in section 5 pointing out a scries of aspects which remain to be 
clarified about this subject. The results are presented in natural units % = c = 1. 



2 The scattering amplitude 

We start with the exact solutions of the free Dirac equation in the de Sitter 
spacetime written in [2]. Let us write the de Sitter line element [1], 

ds 2 = dt 2 -e 2 ^dx 2 , (1) 

where to is the expansion factor and ui > 0. Now we know that defining a spinor 
field on curved spacetime requires one to use the tetrad fields [1] e.jx(x) and 
eP- (x) , fixing the local frames and corresponding coframes which are labelled by 
the local indices (i, 0, ... — 0, 1, 2, 3. The form of the line element allows one to 
chose the simple Cartesian gauge with the non-vanishing tetrad components: 

eg = e-"*; eX = S^, (2) 

so that ejx = e v ^e v and have the orthonormalization properties e^ep = rj/xo, 
e^ep = 8^ with respect to the Minkowski metric r\ — diag(l, — 1, — 1, —1). 

Now let us introduce normalized helicity spinors for an arbitrary vector p by 
notation: £\(p), 

apUp) = 2 P XUP), (3) 

with A = ±1/2 and where a are the Pauli matrices and p =\ p |. For witting 
the solutions of Dirac equation on de Sitter spacetime we set: 

k = — , v± = \±ik (4) 

Then the positive frequency modes of momentum p and helicity A that were con- 
structed in [2] using the gamma matrices in Dirac representation (with diagonal 
7 ) are: 

u P At,x) ~ (27r)3/2 ^ Xe -^/2 H w { 2. e -^ )UP) ) e ' ^ 

where (z) is the Hankel function of first kind. 

Since the charge conjugation in a curved background is point independent 
[7], as in Minkowski case, the negative frequency modes can be obtained using 
the charge conjugation, 



U p -x(x) - V PtX {x) = *7 2 7°(t/p,A(x)) T 



(6) 
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Thus we can restrict ourselves to analyze only the positive frequency modes. 
These spinors satisfy the orthonormalization relations [2]: 

J<Px(-g) 1 '*U RX {x)>y U i?iXI {x)= (7) 

J d 3 x{-gf/ 2 V p , x (x)^V pxi {x)^5 xx ,5 3 (p-p') 

J d 3 x(-. 9 ) 1 / 2 ^ A (x) 7 °t/-, A ,(x) = 0, 

where the integration extends on an arbitrary hypersurface t = const and 
( — .9) 1//2 = e 3uJt . They represent a complete system of solutions in the sense 
that [2] 

f d 3 P J2 [Up,x(t,x-)U+ X (t,x>) + V^ x (t,x)V+ x (t,x')] = e- 3ut 6 3 (x - x>) . (8) 

^ A 

In our calculation we need to find the form of the Coulomb potential on de 
Sitter spacetime which depends on metric. Here we can exploit the conformal 
invariance of the Maxwell equations since the de Sitter metric is conformal 
with the Minkowski one. We can write for the Coulomb field in Minkowski 
spacetime. Then we find the corresponding de Sitter potential: 

A b (x) = H e ~ w *> Ai(x)=0, (9) 

where the hated indices indicate label the components in the local Minkowski 
frames. We also observe that Eq. (9) is just the expression from flat space with 
distances dilated/contracted by the factor e - "*. 

Having the above elements in our mind we can proceed to develop the theory 
for the scattering amplitude on de Sitter spacetime which can be reproduced 
from that in the Minkowski space. The necessary requirements for develop the 
scattering theory is the global hyperbolicity of the space and having a complete 
set of solutions of the free equation for incident field and scattered field (Born 
approximation) with the distinction between positive and negative frequencies. 
It is also important to specify that in our analysis both cases are fulfilled. 

Let ipi( x ) and tpf( x ) be the waves freely propagating in the in and out sectors, 
and we assume that they are both of positive frequency. In direct analogy with 
the Minkowski [3], [4] theory we can define the lowest order contribution in the 
scattering amplitude as follows: 

Ai^f = -ie J d 4 x[~g(x)} 1/2 $ f (x)^A^(x)Mx), (10) 

where e is the unit charge of the field. Our intention is to calculate the amplitude 
of Coulomb scattering for the external field (9) and for initial and final states 
of the form 

ij)i{x) = U p - iiXi (x),if) f (x) = U p - f , Xf (x). (11) 
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Now replacing in Eq. (10) the explicit form of the Dirac spinors and the form 
of the Coulomb potential we observe two remarkable properties. First is that 
the dependence of t, x allows us to split the four dimensional integral into a 
pure spatial integral and a temporal one. The second property is that the pure 
spatial integral have exactly the same form as in the Minkowski space. 

(12) 



/ 



d 6 X- 



-Pll 2 



\x\ \f f 

It is clear that the not so simple part is the temporal integral which contains 
the influence of the gravitational field via the expansion parameter, 



-iaZ- 



8ir\p f - p, 

+ S gn(\ f \ l )e-* k / dzzH^{p f z)H^(p t z) 
Jo 

where we pass to a new variable of integration 



z = 



dzzH^(p f z)H^( Pi z) 
(13) 

(14) 



Note that the integration limits in (13) corresponds to t = ±oo , we assume 
that the interaction extends into the past and future. 

For simplification of our notation we introduce the following quantities: 

A k (f,i) = Al{p f pi) + sgn{\ f \i)Al{p f pi), (15) 

with 



A k (PfPi 







dzzH^{Pfz)H^{ Pl z). 



(16) 



Putting all the above notations together we can write for the scattering ampli- 
tude: 

Ai - f = -hw^w^vt)^®), (17) 

with a = e 2 and stand for the initial and final helicity two spinors. 

The evaluation of the integrals is discussed in Appendix, here we give the 
final result in terms of Dirac delta-function, hypergeometric functions, Euler 
Beta functions and unit step function 



Ak(f, i) = S(p f -Pi) + 9(pi 



sgn(XfXi) 



Pi 



El 

Pt 

8(p f -Pi) + 0(pi-Pf)-f-k [ — 



Pi 



Pi 



i-*(p/-ft)-/;- 

0(Pf-Pi)^f- k 



Ei 
Pf 

El 

pi 



We introduce the following notation for simplify our formula: 



fk(x) =HxY 



2 F 1 (\,l-ik;\-ik; X 2 ) 



cosh(irk) 



~i(x) 



l+ik 



2' 2 

-ik-A 



ik) 

+ ik; X 2 ) 



.(18) 



(19) 



cosh(nk) B{-\,^+ik) 
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where x = |~ or an d f-k (x) is obtain when k — > — fc in (19). 

The above formulas is our result and in remaining paper we will explore some 
of their physical consequences. Before that we must state that the argument \ 
in (19) must be considered in interval < x < 1 (note the argument in the 9 
functions) . This is the domain of convergence of the hypergeometric functions, 
because in the limit \ — ► 1 the 2^1 (a, b, c, x 2 ) functions diverge. The above 
observations are crucial when we compare our expression of scattering ampli- 
tude with the well known expression of scattering amplitude from Minkowski 
spacetime. In our case we see that in (19) appear terms that demand pj ^ pi 
and from this we conclude that there exist a nonvanishing probability for a 
scattering process where the law of conservation for total momentum is lost. 



3 Limit cases 

It is clear that the influence of the de Sitter geometry in the amplitude is con- 
tained in the A/-(f,i) factor which was obtained after integration with respect 
to time variable. This quantity encode via the dependence of the k parameter 
the effect of the expansion of the space on the scattering process. 

We start our analysis having in mind that for ui = we obtain from de Sitter 
metric the Minkowski one and this corresponds to fc = 00. We see that in our 
expression just the terms with Dirac delta-functions will not vanish. Also for 
opposite helicities A OQ (f, i) = and for equal helicities we will obtain: 

A oc (f,i) = 28 X/ x i 5(Pf-Pi)- (20) 
The correspondent factor from Minkowski case for equal helicities is: 

AMinkowski(f, i) = 4#(.E/ - Ei), (21) 

where the notations are : Ef,Ei are the final and initial Minkowski energies 
and E 2 = p 2 + to 2 . Now using the properties of distributions and the fact that 
helicity conservation increase in the ultra relativistic limit (p ~ E 3> to) we 
obtain the amplitude ratio deSitter/Minkowski: 

4S(E f -Ei)-2 m 1 ' 

The fact that we don't recover the Minkowski amplitude can be explained as 
follows. As we point out in the previous section the limits of integration in (13) 
for the time variable extends into the past and future t — > ±00. Observing that 
for t — > 00 the integrand vanishes as e~ wt , while for t — » —00 become divergent. 
In this case we can use the asymptotic formulas for Hankel functions 

H^ 2 \z) = (J^J V2 e ±*(*-"r /a-*/*), z 00. (23) 
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Passing to conformal time t c = — z and using (23) the contributions to the 
amplitude in (16) will be : 

i J dt c e^ p f- pi ^, t c ^-oo. (24) 

The integral (24) produces the distributional term, integrating over an infinite 
interval. Observing that in (24) appears only one semi-infinite axis the ex- 
planation for the 1/2 factor is immediate. To summarize in deSitter case the 
distributional terms will be produced by the evolution in the infinite past when 
the space was contracted. This contraction of the space in the infinite past is 
responsible for the 1/2 factor when we pass to Minkowski limit. 

We will make a few comments about nonconservation of total momentum. 
The difference between (21) and the result in the previous section (18) is that 
the last one does not vanish for pj ^ Pi. As we point out in the previous 
section the total momentum is not conserved in the scattering process on de 
Sitter spacetime. We know that in Minkowski case the conservation of the 
momentum modulus is a direct consequence of the conservation of total energy 
which in turns is a consequence of the translational invariance with respect to 
time, which is lost in de Sitter space. 

Now let us take a closer analysis to helicities from which we deduce a in- 
teresting properties. First we observe that the delta terms 5(pf — Pi) in (18), 
will add for equal helicities and vanish for opposite ones. Now if we admit 
that we deal with a process such that modulus of the final momentum is closer 
to the modulus of the initial momentum pf ~ pi the probability of transition 
between different helicities states will be smaller that that between identical he- 
licity states. From that we can conclude that on de Sitter space is manifesting 
a tendency of conservation for hclicity. 

Let us consider now the opposite case k = which corresponds to a very large 
expansion rate of the space compared with mass of particle. The f±k(x) func- 
tions simplify in this case, observing that for null k the second and the third ar- 
gument in hypergeometric functions became the same and using 2-^1(0, b; b; z) = 
(1 — z)~ a 1 we find 

/(fc=o) (x) - " • (25) 

71" 1 - X 

Introducing (25) in (18) we obtain a pair of terms which summed make the unit 
step functions unnecessary. Giving attention to the delta contributions, one 
obtain for identical initial and final helicities: 

A (f,i) = 2S(p f - Pi ) + ^-^— . (26) 

™ Pi ~ Pf 

For opposite helicities the result is identically null, since for k = it is easy to 
see that Af(pf,pi) exactly cancel. 

An important observation is that (26) can be obtained if we put k = in 
Hankel functions obtaining, 

H^ {z) = T i(£j'\^, (27) 
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in which case the integrals in (16) can be easy evaluated. 

Now we know that the Dirac theory is conformally invariant for a null fermion 
mass, and that our space is conformal with the Minkowski one. The case k = 
is equivalent with a vanishing mass m — and the explanation is immediate. 
Also note that the massless limit is consistent with the identically null amplitude 
in the case when helicity is not conserved. 

In the above limit cases we observe that the amplitude will vanish for oppo- 
site helicities and the conclusion is immediate: in bought cases the total angular 
momentum is conserved. As we point out above this tendency is also manifested 
in the general case. 



4 Calculation of cross section 

We know from the scattering theory that the parameter measured by experi- 
ment in scattering processes is the cross section. In our case the situation is 
complicated and we must reconsider our vision about cross section. The first 
observation is that our amplitude of scattering contain two main contributions 
one linear part (terms with delta Dirac function) and one non-linear part. In 
this situation it is natural to define in general the differential cross section as a 
sum of a linear contribution and a non-linear one. Also it is important to in- 
terpret our non-linear amplitude at square as a probability of scatter when the 
modulus of momentum is not conserved. The definition of linear cross section 
for Coulomb scatter can be reproduced as in Minkowski spacetime 

1 1 dPi 

XiX f 

where the factor i transform one sum in mediation, Pi is the linear probability 
of transition and J is the incident flux in unit of time. When modulus of momen- 
tum is conserved just the linear probability of transition will give contribution 
to cross section. 

For the non-linear differential cross section for Coulomb scatter we will use 
the following definition 

da nl = \Y, -j dp nU (29) 

XiX, 

where P n i is the non-linear probability of transition. In the case when modulus 
of momentum is not conserved just the non-linear probability of transition will 
give contribution to cross section. 

For our further calculations we evaluate the incident flux in unit of time 
using the formula of Dirac current in local frames: 

J" = 4U p - i>Xt (x)rU p - ttXi (x). (30) 

The spatial component of Dirac current can be defined in our case as follows: 
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+HM(*e-*)HW(?ie- ut j\. (31) 

UJ CO 

The spatial components of the four vector current of particles is the incident flux 
in unit of time. We see that this is a quantity that depends of time which is a 
unusual and complex problem. In addition another problem appears, we known 
that when we pass from scattering amplitude to cross section an essential role in 
the treatment when we deal with finitely extended wave is played by the energy 
conservation factor 5(Ef — Ei) . This factor assures that the cross section will 
not depend on the particular form of wave function of the incident particles. 
In our case the cross section calculate with our amplitude will depend on the 
specific form of the incident wave which is another unusual situation. All these 
aspects are consequences of the lost invariance with respect to time in de Sitter 
spacetime. 

For calculating the cross section in this situation we will follow the next 
physical picture. According to the formalism of the scattering theory the am- 
plitude obtained here has to be associated with the following picture: the in 
states is freely propagated from the moment t = back to the infinite past, 
then evaluated in presence of the interaction up to infinite future, and again 
freely propagated back to t — where is projected in out states. In this situa- 
tion we can evaluate Eq. (31) at t = 0, thus obtaining one quantity that is not 
dependant of time, 



J nPi 



4w(2tt) s 



^ 2 _ ) (^)^ ) (-) + < ) (-)< ) (-)l ■ (32) 



When we deal with particles with given helicity in processes of scattering we 
must average upon the helicities of incident particles and sum upon the hclicities 
of the emergent particles. In our case the situation is immediate and we obtain: 

^£ (P5)] 2 = 2. (33) 

For calculation of linear cross section we need to evaluate ^ thus obtaining 
dPi {aZf 



dt 8ir 3 \pf-pi 



8(p f - Pi )d 6 p f . (34) 



Replacing(34) and (32) in (28) writing d 3 pf = PfdpfdQ Pf and integrating with 
respect to pf the factor with Dirac delta- function will be eliminated and in 
addition using that pf = Pi = p we will obtain for the linear part: 



d<n (aZ) 2 



dfl 27rp3sin 4 (f) 



H W { P )H W { P } + H ^ { P )H W { P)\ 1 . (35) 



We see that our cross section have the usual angular dependence as the 
cross section in Minkowski case. In the k = limit we will obtain a much 
simple formula: 

'<M - ^ (36) 



dnj k=0 8p2 si n 4 (|)- 
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Now for calculation of cross section in the Minkowski limit (k = oo) we must 

(1 2) 

approximate the Hankcl functions Hi ' (z) when bought arguments v and z 
have large values which is one of the most difficult problems. In our case one 
of the argument is imaginal and the other real and in the theory of cylindrical 
functions are given just the asymptotic formulas when bought are real or imag- 
inal. Thus this problem of taking this limit is not solved, our calculations in 
this direction shows that this is a very difficult task. For this reason we can't 
evolve the cross section in the Minkowski limit. 

In the non-linear case the integration with respect to pj give an expression 
which is to complicated to be given here. For this reason we restrict to give 
an expression for the differential probability of scatter in this case. Formulas 
necessary to solve integrals with respect to p / for the calculation of cross section 
is given in Appendix B. 



dP {aZf 



M k ( + sgn{\i\ f )B k 

Pi J \Pi 



°(Pf ~Pi)^2 
P f 



\I, ( • sgn(\i\ f )B k 



dtpf (37) 



In the above relation we introduce the following notations: 

B k (x) - h (x) /_* (x) + ft (x) f-k (X) , 

M k ( X ) = \f k (x)\ 2 + \f-k(x)\ 2 • (38) 

In the k = limit we can evolve the non-linear contribution to cross section 
making the integral over final momentum and restrict the limit of integration 
to pfmax because if we take the integrals between and oo our integral diverge. 
In fact this restriction to the limits of integration is suggested by the form of 
the scattering amplitude (26) observing that the final momentum can't exceed 
one maxim value. 

d<T nl \ _ 16(aZ) 2 \Pfmax(2pi - Pfmax) - %Pi (Pi - Pfmax) M Pi -£ max )] 



dtt J k=0 API -Pi\ 4 (Pi- Pfmax) 

(39) 

We see that the cross sections calculated with non-linear part of amplitude 
have much complicated expressions. In the general case as we can see from 
(37) the non-linear cross section will depend on the parameter of expansion and 
will have a complicated expression, which is a sum of expressions of the form 
(44) (see Appendix B). For very small values of u> the spacetime appears as flat 
and contributions that contain negative powers of k can be neglected. In this 
approximation contributions to non-linear cross section will come from terms 
which are proportional to e 2lTk and x 2tk - 

5 Conclusions 

In this paper we investigated the Coulomb scattering amplitude for the Dirac 
field in the expanding de Sitter space. We have complete ignored complications 



9 



due to the ambiguity of the particle concept in a curved spacetime. We consid- 
ered the initial and final states of the field as exact solutions of the free Dirac 
equation in de Sitter space in the momentum and helicity basis. We also found 
that the amplitude depends in an essentially way on the parameter k = m/uj. 
For a vanishing k we recover the expected result due to the conformal invariance 
of the theory in the massless case. 

Also we found that in the vanishing rate of expansion the Minkowski ampli- 
tude is not recovered. We explain this as a consequence of the contraction of 
space in the infinite past. 

Now we want to draw attention to some aspects that remaind untouched 
in the paper and which could be helpful for further investigations. A question 
of principle arises in our paper: how to define the cross section when this is 
a quantity which is dependant of the form of incident wave and in addition 
the incident flux is also a dependant of time quantity. Thus when we calculate 
the cross section we evolve the incident flux at t = 0. This is the only way to 
calculate the cross section in our situation and suggest a possible principle when 
we evolve the cross section in a spacetime where the translational symmetry with 
respect to time is lost. Also we found that the cross section can be defined in 
generally as a sum of a linear and a non- linear part. From our point of view the 
next important step is to calculate one scatter amplitude when the projectile 
particle are at finite distances from target, because in this situation we may 
think to a better connection to experiment. It is also important to say that an 
experiment in our situation will depend of the orientation of the axes of local 
tetrad. 

All these difficult problems that arises in our paper are consequences of the 
lost invariance with respect to time in de Sitter spacetime. 

6 Appendix Arlntegrals of Bessel functions 

We will present here the main steps leading to (18). The formulas that we need 
are [5], [6] 

p 2sin(7r^) 

and integrals of the type Weber-Schafheitlin [5] 

dzz~ s J^az)J,{fiz). (41) 

These integrals can be solved in the assumptions, 

(3 > a > 0, Re(s) > -1, Re(fi + v-s + l)>0 
a > [3 > 0, Re(s) > -1, Re(p + v - s + 1) > 0, (42) 



o 
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bought cases are fulfilled in our analysis. For obtaining (18) we must introduce 
(40) in (16) and using Weber-Schafhcitlin integrals we will arrive at the desired 
result. Now one sees that our integrals in (16) demand s = — 1. This is prob- 
lematic because the integral becomes oscillatory for z — > oo. Then a simple way 
to solve this problem is to consider an parameter s of the form: 

s = -l + e, (43) 

and let in the end e — > 0. No notably differences appear when evaluating the 
integrals in (16) directly with s = — 1 and with [i,v arguments following from 

Let us take a look to the distributional contribution 8{pf — pi). This orig- 
inates in the term 8 {ex. — (3) which is generally present in the integrals (41) for 
s = — 1 (the case of real arguments /i = v is a familiar one). In our case we must 
obtain this term for arbitrary ji, v arguments, because in our case they depend 
of the parameter k. In watt follows we shall briefly indicate below a method to 
obtain it. We start with s = — 1 + e parameter like in (43) and integrate in (41) 
with respect to (3 over the interval a — e < j3 < a + e. Then the coefficient of 
8{cx — [3) will be obtained then as the quantity which is left in the limit e — > 0. 



7 Appendix B:Integrals with hyper geometric func- 
tions 

For calculation of non-linear cross section we solve integrals of the form: 

/*°° Xp 1 7Y^ 

/ HPt - PfhFi{b, c; d; -^) 2 F 1 {e, /; g; -^)p n f dp f 
Jo Pi Pi 

1/ l\-"/ 2 „ [T(-b + c)T(-l + b + e-%)T{g)T(-l + b + f-%) 



2 V P 



-2 ) ptnd) 



T{c)T{-b + d)T{-l + b + g- f )r(e)r(/) 

71 71 71 

xT{l - b + -UF 3 {b, l + b-d,-l + b + e--,-l + b + f--;l + b-c 



n , n 



, , 1 r(-l+6 + /-f)r(l(l-26 + n)) 

6 - 2 ' 1 + h + 9 - 2 ; 1} - V ~pi Pi mn-b + d)r ( -i + b+g -v 

X f£)f(7) 4F 8 (M + 6-d > -- + 6 + e-- ) 

-\+b + f-^l + b-c,\+b-^,- l -+b + g-^l) + T{b-c) 

r(-i + c + e - f )r(. 9 )r(-i + c + / - | )r(i - c+ f) 

r(6)r(- c + d)r(-i + c + 5 -f)r(e)r(/) *wi + c < 

71 71 71 71 

-l + c+e- -,-l + c + /- -;i-6 + c ,c - -,-l + c + 5 - -;1) 



i r(6 - c)r(-i + c + e - !)r(. 9 )r(-i + c + / - §)r(±(i - 2c + n)) 



pf l r{b)r(-c + rf)r(-i + c + . 9 - §)r(e)r(/ ) 



11 



i n \ ft ft 

x 4 F 3 (c, i + c -d,--+c+e--,-- + c+f--;l-b + c,-+c--, 



i /_1 r(-i + 6-f)r(-i + c-f)r(i + f) 

2 +c + ff 2' lj y ^ r(-i + rf-f)r(c)r(6) 

_ , , 1 + n 3 , n 3 , n 3 n, 
x^fe,/,— ,--<*+-;#,-- 6+ -,--c+-;l) 

r(-i + 6-f)r(-i + c -f)r(i + f) F(ef ,n n 
+ r(-i + rf-f)r(c)r(6) 4 F 3 ( e ,/,i + -,2 - d + -, 

5,2-6+-, 2-c+-;l) 

2 

and integrals of the form J °° 6(pi —pf)[ 2 Fi(b, c; d; ^i)] 2 pf 2lk+n dpf which is the 
same type as the above and have a closer result. The above integrals are valid 
in the assumption n > 0. 
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